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Abstract 

Heavy quark production in multijet events at e + e~ colliders is studied at tree 
level. Total production rates are given and compared with the corresponding results 
for massless quarks. A new method of computing helicity amplitudes is briefly 
sketched. 
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Introduction 



The great number of hadronic decays of the Z° observed at LEP provides the opportunity to 
test our understanding of strong interactions in unprecedented detail. Large samples of multi- 
jet events have been accumulated and analyzed Q. Recent advances in 6-tagging techniques 
based on the introduction of vertex detectors and on a refinement of the selection procedures, 
with their large efficiencies and the resulting high purities, have paved the way to the study of 
heavy-quark production in association with light-quark and gluon jets. 

It has been pointed out in series of papers || |], || that the effects of the 6-quark mass are 
substantial and increase with the number of jets. We have studied jet production at tree level 
taking full account of 7, Z interference and of quark masses. These latter reduce the available 
phase space and strongly decrease the emission of gluons collinear with the quark direction, the 
so called "dead cone effect". In particular we have computed: 

• e + e~ — > QQg • e + e~ — > QQgg • e + e" — > Q1Q1Q2Q2 

• e + e~ — ► QQ'j • e + e~ — > QQ^g • e + e~ — > QQjj 

• e + e~ -> QQggg • e+e" -> Q1Q1Q2Q29 

We have used M z = 91.1 GeV, Y z = 2.5 GeV, sm 2 (9 w ) = -23, m b = 5. GeV, a em = 1/128 and 
a s = .115 in the numerical part of our work. 

When the number of Feynman diagrams becomes large it is convenient to compute the 
amplitudes using helicity methods instead of computing directly the amplitude squared. We 
have used two of the most popular formalism || 0] in our calculations. Both methods can 
be easily implemented in a small set of nested subroutines. This however results in computer 
programs which are too slow. For the five-jet case we have resorted to the symbolic package 
Mathematica || to write down the Fortran expression for each helicity amplitude. With this 
procedure we have produced a rather large piece of code, which however runs quite fast, and 
therefeore can be used in high statistics Montecarlo runs. As an example the program for qqggg 
production is about 24,000 lines long, but requires only about 5 x 10~ 2 seconds to evaluate on 
a Vaxstation 4000/90. This is still acceptable but clearly indicates that faster methods are 
needed, as the one we have recently developed and used for computing e + e~ — > bbW + W~ 
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A new method for helicity amplitude calculations 

There are normally three possible ways of evaluating a spinor line. The first consists in reducing 
the expression to a trace The second amounts to writing explicitely, for example in the 
helicity representation 0, the components of the spinors, of the vertex matrices ^'s and of 
the fa and then proceed to the multiplication of the matrices and spinors. The other way 
12, 0] consists in decomposing every j>{ in sums of external momenta ^ and use the relation 
# = Ea U(k, X)U(k, A) + M (with M = +m if U = u, M = -m if U = v) in order to reduce 
everything to the computation of expressions of the type U(ki, \i)i$U(kj, Xj). 



We get a remarkable simplification with respect to the procedures sketched above inserting 
just before every (y^ + /jj) in a spinor line, completeness relations formed with eigenvectors of 
j>i. To do this we must construct spinors U(p, A) which are defined also for p spacelike. With 
this method, in addition to reducing ourselves to the computation of expressions of the type 
Ufa, \i)ijfiU(pj, Xj), we avoid the proliferation of terms due to the decomposition of the in 
terms of external momenta. 

One can easily costruct an example of spinors defined for any value of p 2 and satisfying Dirac 
equation and completeness relation, with a straightforward generalization of those introduced 
in ref.0. One first defines spinors w(k , A) for an auxiliary massless vector k satisfying 

I _|_ \ry 5 

w(ko, X)w(ko, A) = — - — (1) 

and with their relative phase fixed by 

w(k Q ,\) = )Jf 1 w(k Q ,-\), (2) 

with ki a second auxiliary vector such that k\ — — 1, ko-ki — 0. Spinors for a four momentum 
p, with m? = p 2 are then obtained as: 

U (P> A ) = -jt=^= w ( ko > ~ A ) v ivA) = ^j===w{k ,-\) (3) 

u(p, A) = w(k , -A) i±5= «(P, A) = w(k , -A) 4^5= ( 4 ) 

If p is spacelike, one of the two determination of ^/p 2 " has to be chosen for m in the above 
formulae, but physical results will not depend on this choice. 

One can readily check that with the previous definitions, Dirac equations 

^■u(p) = +mu(p) pt>(p) = —mv{p) (5) 

u{p)$ = +mu(p) v{p)f = —mv(p) (6) 

and the completeness relation 

= u(p, \)u(p, A) - v(p, X)v(p, A) 

are satisfied also when p 2 < and m is imaginary. 

Let us now consider the case in which there are only two insertions in a spinor line: 

T < - 2 \p 1 ;r] 1 ;p 2 ;r] 2 ;p 3 ) = U(p u Ai)^i(ff 2 + lMi)$-iU{pzi A 3 ). (8) 

One can insert in eq. (H), on the left of {j> 2 + ^2) the relation (0) and make use of Dirac 
equations to get: 



T (2) = \u{ Vl , X 1 )H 1 u(p 2 , A 2 ) x w(p 2 , A 2 )^ 2 C/(p3, A 3 ) x f 1 + ^) + 

Z \ 777-2 / 

^I7(pi, Ai)^iu(p2, A 2 ) x v(p 2 , X 2 )i 2 U (p 3 , A 3 ) x f 1 - —) (9) 

2 V 777 2 / 



This example can be generalized to any number of insertions and shows that the factors (p'i + Hi) 
can be eliminated, reducing all fermion lines essentially to products of T functions: 

T XlX2 {p 1 ;r ] ;p 2 ) = U(p 1 ,\ 1 )iU(p 2 ,\ 2 ) (10) 

defined for any value of p\ and p\. 

The T functions ( llTf ) have a simple dependence on m\ and m 2 : 

T\ 1 \ 2 (pi,V,P2) = yjprk yjp2-k T Xl x 2 (pi;rj;p 2 ) = (11) 

^AxAa(Pi;»7;P2) + MiB Xl \ 2 (pi,V,P2) + M 2 C MX2 (p 1 ;r);p 2 ) + M 1 M 2 D Xl)vi (p 1 ;ri;p i ) 

where 

Mi = +rrii if U(p h A») = A») M { = -mi if Z7(p i; A») = u(p i; A»). (12) 

The functions A, 5, C, -D are independent of mi and m 2 and of the u or v nature of U(p\, Ai) 
and U(p 2 , X 2 )- Every T^ n \ for any number of insertions turns out to be of the form (|T0|) and 
this greatly simplyfies the rules for evaluating spinor lines || . 



Results 

A selection of our results is shown in fig.l through 3, and we refer to the original papers for 
more details. In fig. lb we present the cross sections for e + e~ — > qqg and e + e~ — > qqgg with 
q = d, b as a function of y cnt for both the JADE [I3| and DURHAM [H| definitions of y at 



LEP I. The ratio of massive to masslees cross section is shown in fig. la. For small y® nt the cross 
section for bbg is almost 20% smaller than for ddg. As expected the ratio becomes closer to one 
for larger y cnt , but for y cut as large as .2, still R^ d = a (bbg) /a (ddg) < .96 in both schemes. 

Jet-shape variables have been extensively studied as a tool to determine a s and as a testing 
ground for the agreement between data and the standard description of strong interactions. In 
the ranges used for measuring a s , the ratio of massive to massless tree-level predictions can 
significantly differ from unity and it depends both from the variable and from its actual value. 
We have compared at 0(a s ) the ratio Ryz, which is obtained from the full matrix element, 
with Ry, the ratio which results neglecting the Z° (as in JETSET), for Thrust, Oblateness, C- 
parameter, Mh and Mb- The difference between R^z and i? 7 turns out to be about 1.2 x 10~ 2 , 
almost independent of the particular variable and of its specific value. As an example, in fig. 2 
we show both R jZ and _R 7 for the M H , M D , C and O distributions. 

In fig. 3 we show the ratios cr(2b3g)/cr(2d3g) (continuous line), a(2u2bg) / a(2u2dg) (dashed 
line) and a(2d2bg)/a(2d2sg) (dotted line) in the two recombination schemes as a function of 
y cut . These curves confirm our previous conclusions that mass effects increase with the number 
of final state light partons. The ratio for the dominant 2g3g production process is equal to .58 
at y cu t = -001 in the DURHAM scheme and to .67 at y cu t = .005 in the JADE scheme. It is 
even smaller for the processes with four quark jets in the final state. This corresponds to a 
6-i-8% decrease in the predictions for the total five-jet cross section. 

When our results for jet production are compared with the data it is to be remebered that 
we have used a s = .115 which corresponds to Q 2 = M 2 Q with A-^ = 200 MeV with five active 



flavours. The analysis of shape variables and jet rates to 0{a 2 s ) has shown that, in order to get 
agreement between the data and the theoretical predictions, a rather small scale for the strong 
coupling constant has to be chosen ||. 

The flavour independence of the strong coupling constant has been investigated by several 
groups |L5| using both jet-rates and shape variables. Mass corrections have played an important 
role in this study. We look forward to more detailed analyses of QCD in the heavy quark sector. 
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Figure Captions 

Fig. 1 In the upper part we show Rf 1 = a(ccg)/a(uug), = a(bbg) / 'a(ddg) and R b 4 d = 
a(bbgg)/cr(ddgg) as a function of y cut for y J (dashed) and y D (continuous). In the lower 
part we present the cross sections for e + e~ — > bbg (continuous), e + e~ — > ddg (dashed), 
e + e~ — > bbgg (chain-dotted) and e + e~ — > rfJg^ (dotted) as a function of y cut for both 
definitions of y at y/s = 91.1 GeV. 

Fig. 2 The ratios R = da(bbg)/dM/da(ddg)/dM for M = M H /y/s and M = M D /y/s (lower 
part) and for M = C —parameter and M = Oblateness (upper part) from the full matrix 
element (continuous) and from the photon contribution alone (dashed) at y/s = 91.1 GeV. 

Fig. 3 Ratio of massive to massless cross sections for a(2b3g)/a(2d3g) (continuous line), 
a(2u2bg)/a(2u2dg) (dashed line) and a{2d2bg) / a{2d2sg) (dotted line) in the JADE and 
DURHAM recombination schemes as a function of y cut at \fs = 91.1 GeV. 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405384vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405384vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405384vl 
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